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J*V^ Abstract. We describe the first results from our project aimed at large-scale calculations of 

jj^ i NLTE abundance corrections for important astrophysical atoms and ions. In this paper, the 

focus is on Fe which is a proxy of stellar metallicity and is commonly used to derive effective 
temperature and gravity. We present a small grid of NLTE abundance corrections for Fe I lines 
and discuss how NLTE effects influence determination of effective temperature, surface gravity, 
O ' and metallicity for late-type stars. 

CO 

1. Introduction 

Iron is a key element in stellar astrophysics. The very complex atomic structure of its lowest 
ionization stages, Fe I and Fe II, gives rise to a wealth of spectral lines all across the spectrum 
of a typical late-type star. This atomic property, coupled to a relatively large abundance makes 
Fe a reference for spectroscopic estimates of stellar parameters using the method of excitation- 
Q\ ■ ionization balance. 

The classical implementation of this method in spectrum analysis codes involves three 
assumptions: local thermodynamic equilibrium (LTE), hydrostatic equilibrium, and ID 
geometry. These approximations strongly reduce the complexity of the problem, thus permitting 
analysis of very large stellar samples in short timescales. Yet, in the conditions when the 
breakdown of ID static LTE models occurs the inferred stellar parameters suffer from large 
systematic biases. To assess the latter, more physically realistic modeling is necessary. 

Studies of NLTE effects on the Fe I/Fe II level populations for FGK stars trace back to Athay 
&: Lites (1972). Since then, vast amount of work has been performed in this field demonstrating 
that NLTE effects in the excitation-ionization balance of Fe I/Fe II are significant and can not be 
ignored even in the analysis of solar- type stars (Mashonkina et al. 2011 and references therein). 
Yet, despite major efforts aimed at understanding how non-equilibrium thermodynamics affects 
the line formation of Fe, there has never been an attempt to accurately quantify these deviations 
in a systematic manner across a wide range of stellar parameters, and to apply them to a large 
stellar sample. 

Here we present a new NLTE model atom of Fe I/Fe II to be applied in large-scale calculations 
of NLTE abundance corrections for late-type stars. We discuss the NLTE effects influencing 
atomic level populations at the typical conditions in their atmospheres and provide a small 
grid of NLTE corrections for the Fe I lines. The model has been tested on a number of 



well-studied stars with independently-determined parameters, including metal-poor giants and 
turnoff stars. These tests performed with classical ID hydro-static model atmospheres and 
averages of 3D hydro dynamical simulations of stellar convection will be presented in Bergemann 
et al. (in preparation). A complete grid of NLTE abundance corrections computed with 
multi2.3 statistical equilibrium code will be presented in Lind et al. (in preparation). 

2. Methods 

2.1. Model atmospheres and codes 

The calculations presented here were performed with ID LTE plane-parallel mafags-ODF models 
(Fuhrmann et al. 1997, Grupp 2004). In these models, convective energy transport is accounted 
for using the mixing-length theory of Bohm-Vitense (1958) with the mixing length parameter 
Omit set to 0.5. Line opacity is represented by all elements up to Ni, and various diatomic 
molecules (H2, CH, CO, TiO, etc). This accounts for nearly 20 million atomic and molecular 
lines. The reference solar abundances were compiled from various literature sources, giving 
preference to the determinations by the Munich group. The model atmosphere provide partition 
functions and partial pressures of all relevant atoms and molecules, which are then used in 
detailed line formation calculations. 

NLTE statistical equilibrium was computed with a revised version of the detail code (Butler 
& Giddings 1985), which solves multi-level NLTE radiative transfer with a given static ID model 
atmosphere. The last version of the code is based on the method of Accelerated Lambda Iteration 
(ALI). LTE and NLTE line formation calculations with departure coefficients from detail were 
performed with the revised version of the spectrum synthesis code SIU (Reetz 1999). The major 
update in the code relates to the automated computation of NLTE abundance correctional, 
which generally proceeds by constructing the grids of LTE and NLTE line equivalent widths for 
a range of stellar parameters and interpolating in the corresponding curves-of-growth for the 
input observed equivalent widths or abundances. 

2.2. Model atom 

The NLTE model atom of Fe was constructed using all available atomic data from the Kuruc^] 
database, which also includes laboratory data from the NIST compilation^. All predicted 
energy levels of Fe I with the same parity above -Blow > 5.1 eV were grouped into super- 
levels. Transitions between the components of the super-levels were also grouped. The total 
transition probability of a super-line is a weighted average of logg/'s of individual transitions. 
Thus, in the final atomic model, the number of energy levels is 296 for Fe I and 112 for Fe II, 
with uppermost excited levels located at 0.03 eV and 2.72 eV below the respective ionization 
limits, 7.9 eV and 16.19 eV. The model is closed by the Fe III ground state. The total number 
of radiatively-allowed transitions is 16 207 (13 888 Fe I and 2 316 Fe II). 

Photoionization cross-sections for 136 states of Fe I were taken from Bautista (1997) and 
the hydrogenic approximation was adopted for all other levels. The rates of transitions induced 
by inelastic collisions with free electrons (e~) and H I atoms were computed using different 
recipes. For the states coupled by permitted b-b and b-f transitions, we used the formula of 
van Regemorter (1962), respectively, Seaton (1962) for e~. To derive the rates of transitions 
due to inelastic collisions with H I atoms, we used the Drawins formula (Drawin 1969) in the 
version of Steenbock & Holweger (1984). Also, the states are connected by forbidden transitions 
induced by inelastic collisions with e~ and H I atoms. These are computed using the formulae 

1 NLTE abundance correction is denned to be the difference in abundance required to match NLTE and LTE 
line profiles or equivalent widths 

2 http: / /kurucz. harvard.edu/atoms. html 

3 http://www.nist.gov/pml/data/asd.cfm 
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Figure 1. Left panel: Grotrian diagram of the Fe I atom. Right panel: comparison of quantum- 
mechanical rate coefficients for e~ -induced forbidden transitions between the lowest energy levels 
of Fe II with the classical Allen's (1973) recipe. 

of Allen (1973) and Takeda (1994). Quantum-mechanical calculations exist only for e~ -induced 
forbidden transitions between the 16 lowest Fe II energy levels (Ramsbottom et al. 2007). These 
data are compared with the Allen's (1973) recipe in Fig. [TJ 

The influence of e~ and H I collisions on the statistical equilibrium of Fe was carefully 
investigated by performing test calculations with various scaling factors to the above-mentioned 
formulae. The NLTE synthetic profiles were furthermore compared with the observed stellar 
spectra to check how well the test model atoms satisfy the constraint of ionization-excitation 
balance of Fe I/Fe II under restriction of different stellar parameters. The final choice of scaling 
factors will be discussed in Bergemann et al. (in preparation). 

3. Statistical equilibrium of Fe 

3.1. NLTE effects 

Fig. [2] shows level departure coefficients, 6jJl as a function of continuum optical depth logrsooo 
for a number of model atmospheres computed with different stellar parameters. The latter 
are representative of stars typically used in Galactic chemical evolution studies. Only selected 
energy levels of Fe I and Fe II typical for their depth-dependence are included in the plot: a 5 D 
(ground state of Fe I), z 7 D° (2.4 eV), e 5 D (5.4 eV), t 3 R a , a 6 D (ground state of Fe II), and 
z 4 D° (5.5 eV). All these levels give rise to radiatively-permitted transitions, which are used in 
the detailed abundance analysis of the Sun and a number of metal-type stars. 

The common result for all the models shown in Fig. [2] is that in the optically thin atmospheric 
layers the majority of Fe I levels are underpopulated compared to LTE, bi < 1. The Fe II number 
densities remain close to LTE values, bi ~ 1, and a minor overpopulation of the Fe II ground 
state develops only close to the outer atmospheric boundary. 

Deviations from LTE in the distribution of atomic level populations arise because mean 
radiation field, J u , at different depths and frequencies is not equal to the Planck function, 
B u [T c (t)], which is defined by the local temperature T e (r) at each depth. For Fe I, excess of 
the mean intensity over Planck function in the UV continua leads to over-ionization, which sets 
in as soon as the optical depth in the photoionization continua of low-excitation Fe I levels 
falls below unity, i.e. logrsooo ~ 0. In the outer layers, ionization balance is dominated by the 



4 Departure coefficient is denned to be the ratio of NLTE to LTE atomic level populations, bi = nf LTE /ni' 
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Figure 2. Departure coefficients of selected Fe I and Fe II levels for a number of model 
atmospheres from the grid; stellar parameters are indicated in each figure. 



Fe I levels with excitation energies at 4 — 5 eV. In the infra-red continuum, J v < B v leads to 
over-recombination, which is very efficient for our atomic model with only 0.03 eV energy gap 
of the upper Fe I levels from the Fe II ground state. 

Due to an extremely tight radiative and collisional coupling of atomic energy levels, excitation 
balance of Fe I is also affected by radiative imbalances in line transitions. These include radiative 
pumping by super-thermal UV radiation field of non-local origin at the line frequencies, as 
well as photon suction driven by photon losses in large-probability transitions between highly- 
excited levels. These processes leave a characteristic imprint on the behavior of 6j-factors in the 
outer atmospheric layers. In Fig. this is seen as the onset of deviations from the relative 
thermal equilibrium between the Fe I levels, which occurs, depending on the model metallicity, 
at logT5ooo ~ — 1... — 2. The NLTE effects on the line formation are different for the strong and 
weak Fe I lines. The former are mostly shaped by deviations of their source function from the 
Planck function, which makes a clear effect in the line cores. Weak lines are pre-dominantly 
controlled by their opacity, thus the major change is due to the shift of their optical depth scale 
in NLTE. 



3.2. NLTE abundance corrections 

We computed NLTE abundance corrections, hereafter Anlte> to Fe I lines for a small grid of 
model atmospheres covering metallicity range —3.6 < [Fe/H] < +0.3. The results are presented 
in Fig. for the parameters representative of giants (left panel) and dwarfs (right panel). 
Our sample of Fe I lines includes transitions between to 5 eV levels, and the line equivalent 
widths range from 5 mA to few A depending on stellar parameters. In the plot, each symbol 
corresponds to a mean Anlte obtained by averaging individual NLTE corrections for all Fe 
I lines in the sample, and the errors bars correspond to the line-to-line variations in A for 
each model atmosphere. The former quantity roughly demonstrates how much the Fe I/Fe II 
ionization balance, and, thus, determination of surface gravity, is affected by NLTE, whereas 
the latter helps to assess the influence of NLTE on the excitation balance and determination of 
effective temperature. 
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Figure 3. NLTE abundance corrections for different stellar parameters. Right panel: all model 
atmospheres adopt log g = 4.6. 

NLTE corrections increase with decreasing [Fe/H] and log g, thus ionization balance achieved 
assuming LTE leads to progressively underestimated gravities and metallicities. In particular, 
at T e ff = 5000 K (Fig. [3l left panel), the mean NLTE corrections are twice as large in the 
log g = 2.2 model compared to the log g = 3.6 model, reaching +0.4 dex at [Fe/H] ~ —3.5. 
It is also interesting that NLTE corrections for dwarf models critically depend on T e g. Fig. 



[3] (right panel) shows that for the coolest models, T g ~ 5000 K, one may safely adopt LTE 
approximation down to the lowest metallicities. However, LTE leads to severe errors in gravity 
estimates for warm turn-off stars, T > 6000 K. 

On the other side, non-negligible error bars indicate that there are line-to-line variations of 
NLTE abundance corrections for the same model atmosphere. The spread of Anlte among 
Fe I lines with different equivalent widths W\ and excitation potentials E\ ow strongly increases 
with decreasing metallicity for giants, and less so for dwarfs. Furthermore, this scatter is not 
random, but corresponds to a trend of Anlte with E\ ow , as shown in Fig. HI The figure 
demonstrates variation of NLTE corrections with lower level excitation potential for the models 
with parameters: T eff = 5000, log g = 2.2, [Fe/H] = -2.4 (left panel) and T eff = 5800, log g = 4.6, 
[Fe/H] = —2.4 (right panel). Stronger lines, W\ > 60 mA, tend to experience smaller departures 
from LTE compared to weaker lines, but they also demonstrate systematically increasing Anlte 
with increasing E\ ow . No such effect is seen for the weak lines, W\ < 60 mA, which, for some 
exceptions, have very similar NLTE corrections independent of their excitation potential. 

Thus, our results suggest that LTE T e g determinations become very inaccurate at [Fe/H] 
< — 2 and there are systematic effects in the T e g- scale of dwarfs and giants. 
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Figure 4. NLTE abundance corrections for the Fe I lines as a function of lower level excitation 
potential. Red and blue colors separate the spectral lines with W\ greater and less than 60 mA, 
respectively. Stellar parameters are indicated in the plots. 



4. Conclusions 

We have constructed a new NLTE Fe model atom using the most up-to-date theoretical and 
experimental atomic data. The model has been tested on a number of well-studied late-type stars 
with parameters determined by other independent methods (Bergemann et al, in preparation). 

Our preliminary results for the NLTE effects on Fe I and Fe II in the atmospheres of late- 
type stars are qualitatively consistent with other studies. Kinetic equilibrium of Fe, which is a 
minority ion at the typical conditions of these cool and dense atmospheres, favors lower number 
densities of Fe I compared to LTE. The number densities of Fe II are hardly affected by NLTE. In 
general, this leads to a weakening of Fe I lines compared to LTE, which, in turn, requires larger 
Fe abundance to fit a given observed spectral line. The magnitude of departures and NLTE 
corrections critically depends on stellar parameters. Whereas NLTE abundance corrections for 
the Fe I lines do not exceed ~ 0.05 dex for solar metallicity dwarfs, they can be as large as 0.6 
dex for metal-poor giants. 

We find that both excitation and ionization balance is affected by NLTE effects on Fe I. 
NLTE corrections to Fe I lines increase with decreasing [Fe/H] and log g, thus ionization balance 



achieved assuming LTE leads to progressively underestimated gravities and metallicities. On 
the other side, Fe I levels are also not in thermal equilibrium with respect to each other, and 
for the stronger lines NLTE abundance corrections show trends with line excitation potential. 
In particular, we find that LTE T c g determinations become very inaccurate at low metallicities 
and significant systematical effects in the T e g scale of dwarfs and giants can be expected. 

The NLTE abundance corrections will be publicly available through an interactive online 
database, which is currently under construction. They can be used to retrieve NLTE abundances 
and abundance corrections for input equivalent widths. We also consider an option to provide 
functional dependency of NLTE corrections on line parameters for a grid of stellar parameters 
(T e ff, logg, [Fe/H]). This is a more convenient solution for implementation in the automated 
codes used in large-scale spectroscopic studies, such as follow-up spectroscopy for Gaia targets. 
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